Density, short-range order and the quark-gluon plasma 
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We study the thermal part ol the energy density spatial correlator in the quark-gluon plasma. We 
describe its qualitative form at high temperatures. We then calculate it out to distances » 1.5/T 
in SU(3) gauge theory lattice simulations for the range of temperatures 0.9 < T/T c < 2.2. The 
vacuum-subtracted correlator exhibits non-monotonic behavior, and is almost conformal by 2T C . Its 
broad maximum at r ~ 0.6/T suggests a dense medium with only weak short-range order, similar 
to a non-relativistic fluid near the liquid-gas phase transition, where n/s is minimal. 

PACS numbers: 12.38.Gc, 12.38.Mh, 25.75.-q 



Hydrodynamics calculations [l[ successfully described 
the pattern of produced particles in heavy ion collisions 
at RHIC [2]. This early agreement between ideal hy- 
drodynamics and experiment has been refined in recent 
times. On the theory side, the dissipative effects of shear 
viscosity 77 have been included in full 3d hydrodynam- 
ics calculations p| 0, 111 and the sensitivity to initial 
conditions quantitatively estimated [6J for the first time. 
On the experimental side, the elliptic flow observable v?, 
which is sensitive to the value of r\ in units of entropy den- 
sity s, is now corrected for non-medium-generated two- 
particle correlations Q- The conclusion that r\js must 
be much smaller than unity has so far withstood these 
refinements of heavy- ion phenomenology Q . 

The smallness of 77/s was turned into the statement 
that the quark-gluon plasma (QGP) formed at RHIC is 
the "most perfect liquid known in nature" Q. A gen- 
eral question then comes to mind: what observable can 
be used to characterize the liquid nature of a system de- 
scribed by a quantum field theory 9]? And secondly, 
what is the QCD prediction for that observable? This 
leads us to remind ourselves what the defining property of 
an ordinary liquid is. Surely the everyday-life notion that 
a liquid "has a definite volume, but no definite shape" is 
inadequate in the present context. 

The two-body density distribution p(ri,r2) = g( r )p 2 
of an ordinary substance (such as water) of density p 
behaves qualitatively differently in the solid, liquid and 
gas phase (see for instance [lOjj ) . The radial distribution 
function g(r) characterizes the average density of parti- 
cles at distance r from an arbitrarily chosen particle. In 
a dilute gas, g(r) is essentially equal to 1 for r greater 
than the size of a molecule. In a liquid on the other hand, 
g(r) vanishes at small r, a reflexion of the short-distance 
repulsion between molecules. The function then rises and 
typically exhibits several gradually damped oscillations 
around unity. This reflects the "short-range order" in 
the fluid, namely the coherent motion of closely packed 



molecules up to distances a few times the molecule size. 
Over longer distances, this ordering is lost. Only a per- 
fect crystal at low temperatures exihibits truly long-range 
order. 

In quantum field theory, particle number is not (nec- 
essarily) conserved, so it is not immediately clear which 
spatial correlator is the closest analogue of the two-body 
density distribution in non-relativistic systems. In QCD, 
the only conserved quantities are energy, momentum, the 
quark numbers and the (non-singlet) axial charges. The 
energy density is an order parameter in the theoretical 
limit of a large number of colors N c , which puts it in natu- 
ral correspondence with the density p of a non-relativistic 
fluid. Hence our strategy to study the spatial correlator 
of the energy density. 

A peculiarity of quantum field theory is that energy 
density correlations are present even when the average 
energy density is zero, i.e. at T = when the partition 
function is saturated by the quantum vacuum. This cor- 
relation has to be strong at short distance in an asymp- 
totically free theory, (T O o(0)T O o(r)) ~ r~ 8 , on dimen- 
sional grounds. Based on the Kallen-Lehmann represen- 
tation, it is monotonically decreasing in r at any T. 

Therefore, in order to isolate the thermal effects, we 
shall consider the subtracted correlator 

G ee (T,r) = <T oo (0)Too(0,r)) T - (r oo (0)T O o(0, r)) T=0 , 

(1) 

In the Euclidean SU(iV c ) gauge theory, the energy den- 
sity operator, 
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is composed of two terms which are separately scale- 
independent operators. Here e is the energy density, p 
the pressure and (3(g) — —bg 3 + . .. the beta function, 
with b = sf^yi ■ The expectation value of the 'naive' 
energy operator #00 is proportional to the entropy den- 
sity, while the expectation value of the trace anomaly 
directly measures the deviation from the conformal limit 



2 



where e = 3p (we implicitly add a constant to Too such 
that (Tqo) vanishes in the vacuum). We can thus investi- 
gate separately the G ss and Gee correlators, defined by 
replacing Too respectively by |0qo and in Eq. [TJ 



High-temperature behavior 

The one-loop expression for the two correlators in D 
space-time dimensions is 



(0oo(O)0oo(aO)iL = (8irba s )- 2 (6(Q)0(x)) 
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where cLa — N% — 1 and a;r„i = (y + xo,x), while 
(6*oo(0)6*(0, r))iL vanishes identically. The to = n = 
term gives the zero-temperature expression. 

Let us consider the high-temperature regime, T 3> 
g 2 (T)T ^> T c . In that case, the behavior of G ee is char- 
acterized by three regimes: 

Tr <C g~ 2 (T): both the zero and high temperature cor- 
relators are well described by the one-loop formula, 
up to small radiative corrections. For r — > 0, 



G ee (T,r) ~ 
(see Fig. El). 
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for Tr > 0.568, G ee > e 2 (T) 



.9 2 (T) <Tr< the zero temperature correlator is 

still described by perturbation theory, , but the 
high-T correlator is exponentially screened, so nec- 
essarily G ee < e 2 (T). 

Tr > ^: both the zero and high temperature correlators 
are exponentially screened. In the former case, the 
relevant mass is M4 ~ 5.3T C [Tll.[l3|. corresponding 
to the lightest scalar glueball mass in D = 4, while 
at high temperatures, dimensional reduction takes 
places, and the screening mass is (M^/g 2 ) ■ g 2 (T)T 
with M 3 /g 2 ~ 2.4 [11]. It is therefore clear that 
screening of the energy density is stronger at high 
temperatures, and hence G ee asymptotically ap- 
proaches e 2 (T) from below at a rate e~ MiT . 

In summary, at high temperatures G ee (T, r) — e 2 (T) van- 
ishes at least at two finite distances r: the first time at 
Tr w 0.568, and the second time at Tr = 0{g- 2 {T)). 
To investigate the function G ee (T, r) at temperatures ac- 
cessible in heavy-ion colliders, we perform lattice simula- 
tions in the region 0.9T C < T < 2T C . 



Operator product expansion 

From expression ^ and from known results [TH, we 
can obtain the leading terms in the operator-product ex- 
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FIG. 1: The correlators at T — at different lattice spacings 
(ro ~ 0.5fm [l||; the line is to guide the eye). 
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FIG. 2: Gee (Eq. [lj across the deconfining phase transition, 
at j3 = 6.018. At T = 1.24T C we check for discretization errors 
by also showing data from /3 = 6.200 (filled squares). 



pansion (OPE) of these correlators, 
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where a; = (0, r) and r~ 2 terms and softer have been 
omitted. The Wilson coefficients of the operators 1, #00 
and are a least of 0(a 2 ) for the product OO (O)0(O, r). 
The usefulness of the OPE in this context arises because 
of the exact cancellation of the r~ 8 term in the difference 
between finite T and T = correlators. Using Eq. [4] we 
obtain the short-distance behavior 

G(T,r) ~ -0^ + O( a °) e -^ + O(a s r-\r- 2 ) (5) 
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G ee (T,r)/(d A T 8 ) 
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FIG. 3: G e e at three temperatures, with iV T = 6. The curve 
corresponds to the correlator in the Stefan-Boltzmann limit 
T -> oo. 



3 
2.5 

2 
1.5 

1 h 
0.5 



-0.5 





9G ss /(16d A T 8 ) 



G flfl /(d A T ). • 



0.5 



1 

Tr 



1.5 



FIG. 4: The N T = 6 entropy-entropy (G ss ) and action-action 
(Gee) correlators at 1.24 (circles), 1.65 (squares) and 2.20T C 
(triangles). 



Numerical results 

We now turn to a calculation of G ee (T, r) in the SU(3) 
gauge theory, i.e. in the plasma of gluons, using lat- 
tice Monte-Carlo techniques on a (1/T) • L 3 lattice. We 
employ the (isotropic) Wilson action [15[ and the 'once- 
HYP-smeared' 'clover' discretization of £?no and 9 devel- 
oped in pjj]. The variance of (#oo) was shown pj} to be 
reduced by almost two orders of magnitude as compared 
to the simplest 'plaquette' discretization. This technical 
improvement allows us to obtain a signal for G ee — e 2 (T) 
out tor - 1.5/T. 

Figure [1] displays the relevant correlators at zero tem- 
perature. They fall off monotonically as r~ 8 at short dis- 
tance and exponentially at large distance. Note that the 
trace anomaly correlator, while 0(a 2 s ) at short distances, 
dominates at large distances. For r/a > 3 comparison of 
the data obtained at three lattice spacings shows that 



discretization errors are under control. 

Figure [5] shows the qualitative change of G ee across 
the deconfming phase transition. At T = 1.24T C the 
functions obtained from these two lattice spacings are 
in qualitative agreement, a non-trivial check, given the 
large cancellation between the finite and zero tempera- 
ture correlators. The function G ee (r) is large and neg- 
ative at short-distances, crosses the asymptotic value 
e 2 (T), reaches a maximum and presumably decreases 
monotonically after that. Although the signal becomes 
too small to tell beyond Tr ~ 1.5, this is plausible in 
view of the small value of the thermal screening mass 
(see next section). 

The 1/r 4 short-distance divergence is not visible be- 
low T c , a fact that the OPE and the smallness of (e,p) 
in the confined phase easily account for. We note that 
G ee /e 2 (T) reaches around r = 0.6/T a maximum which 
is larger at 1.08 than at 1.24T C . We understand this in 
terms of the larger fluctuations present near the (weakly) 
first order phase transition. 

Figure [3] shows the temperature dependence of G ee 
up to 2.2T C . The position of the maximum remains 
r max « 0.6/T, and G ee (r max , T)/e 2 (T) decreases slowly 
as the temperature rises. The curves at 1.65 and 2.20T C 
exhibit near-conformal behavior (i.e., G ee is essentially 
a function of Tr), however the asymptotic approach to 
e 2 (T) has the opposite sign, as a study of screening 
masses shows. 

Figure [5] shows separately the entropy density correla- 
tor and the trace anomaly correlator. The former is qual- 
itatively similar to the energy density correlator, while 
the latter has a rather featureless monotonic behavior. 



Screening masses 

The asymptotic large-r behavior of G ee is dictated by 
the smallest screening mass that Too couples to. This is 
the state invariant under all the symmetries of a constant 
'z-slice' [l9|,[2(j, which has a volume (1/T) xLxL. From 
D = 4 simulations, we obtain directly 



M(T) 
Mi 



0.630(14), 0.906(20), 1.276(32) (6) 



respectively at 1.24 (N T = 8), 1.65, and 2.20T C (N T = 6). 
So it is only at T* = 1.790(36)T C that the thermal screen- 
ing starts to exceed the D — 4 glueball mass. In parti- 
cular, for T > T*, G ee {r,T) approaches its asymptotic 
value from below, and therefore crosses e 2 (T) twice. 



Comparison with non-relativistic systems 

The radial distribution function g(r) of a simple non- 
relativistic liquid, such as 36 Ar at 85K [2l|, exhibits 
several very pronounced peaks above 1. In particular 
g(r) — 1 is of order unity at the first peak. However, it is 
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known [2^ | that 77/s is minimal near the liquid-gas phase 
transition, and becomes large both at low and high tem- 
perature. A highly ordered mesoscopic scale favours the 
transport of momentum, because the holes between the 
closely packed molecules then play the role of quasiparti- 
cles with a long mean free path (an argument attributed 
to Enskog [22]). Heating up the liquid has the effect of 
reducing the amplitude of the peaks in g(r), until they 
disappear completely once the system is in a dilute gas 
phase. Thus the regime where 77/s is minimal is the one 
where g(r) has only few, small oscillations around unity. 

We have found that G ee has exactly one broad peak 
above e 2 (T), exceeding that value by about 5 — 15% (the 
figure decreases slowly with temperature). By analogy 
with non-relativistic fluids, it is tempting to see a relation 
between this fact and the small value obtained for the 
shear viscosity 23] in the same range of temperatures, 
77/s < 1.0. 

Conclusion 

We have calculated non-perturbatively the thermal 
part G ee of the energy-density spatial correlator in the 
plasma of SU(3) gluons, in the range of temperatures 
0.9 < T/T c < 2.2, and found the qualitative high-T 
behavior. It diverges as — r~ 4 near the origin as dic- 
tated by the OPE, and above T c reaches a maximum 
at r max « 0.6/T typically 10% above its asymptotic 



value e 2 (T). For T < 1.8T C it asymptotically approaches 
that value from above, and at higher temperatures it ap- 
proaches from below. While the appearance of G ee is 
rather different from the radial distribution function of a 
typical liquid, we pointed out that it is precisely when the 
short-range order is weak that the ratio 77/s is minimal. 

We note that the radial distribution function of 
monopoles has been computed in SU(2) gauge the- 
ory [24[ , and has a similar shape to G ee . Color charge and 
monopole-antimonopole correlators that look alike have 
also been found in models [25l |. It would be interesting 
to see whether such models can reproduce G ee . 

A straightforward extension of this work is to consider 
the full spatial correlations of the energy momentum ten- 
sor (TpyTpv), in SU(3) gauge theory and in full QCD. As 
a benchmark it would be helpful to know the form of 
these correlators in the strongly coupled J\f = 4 SYM 
theory, which is known to be an excellent fluid from the 
smallness of 77/s [26[ and from its spectral functions [27| . 
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